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ABSTRACT 

In an earlier paper, the effect of the nearby known supernova remnants (SNRs) on the lo- 
cal cosmic -rays (CRs) was studied, considering different possible forms of the particle in- 
jection time. The present work is a continuation of the previous work, but assumes a more 
r*j \ realistic model of CR propagation in the Galaxy. The previous work assumed an unbounded 

Q-r three-dimensional diffusion region, whereas the present one considers a flat cylindrical disc 

bounded in both the radial and vertical directions. The study has found that the effect of the 

■ vertical halo boundary H on the local SNR contribution to the observed CR anisotropy is neg- 
ryj \ ligible as long as H > 2kpc. Considering the values of the halo height H > 2kpc obtained 

by different authors, the present work suggests that the study of the effect of local sources 
on the CR anisotropy can be carried out without having much information on H and hence, 
CN| ■ using the much simpler three-dimentional unbounded solution. Finally, the present work dis- 

cusses about the possibility of explaining the observed anisotropy below the knee by a single 

■ dominant source with properly chosen source parameters, and claims that the source may be 
\ an undetected old SNR with a characteristic age of ~ 1.5 x 10 5 yr located at a distance of 

■ ~ 0.57kpc from the Sun. 
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1 INTRODUCTION 
•rH . 

■ There is a wealth of literature available on the highly isotropic nature of cosmic-rays (CRs) observed at the Earth (see e.g. the references 
\ given in Thoudam 2007, hereafter Paper I). The CR anisotropy amplitude is only ~ (10~ 4 — 10~ 3 ) in the energy range of (10 11 - 10 15 )eV 

- - < (Guillian et al. 2007 and references therein) with the phase (direction) mainly found in the outer Galaxy, particularly in the second quadrant of 
the Galaxy. The possible explanations for the anisotropy are generally beleived to be the global diffusion leakage of CRs from the Galaxy, the 
random nature of the CR sources in space- time and the effect of the local sources. In Paper I, the effect of the known local supernova remnants 
(SNRs) has been studied in detail by giving more emphasis to the particle release time. The study found that the observed anisotropy data 
favour the burst-like injection model if particles are released from the sources at an age of ~ (2 — 5) x W 4 yr. The continuous injection model 
gives an anisotropy which is too large to explain the observed data. However, Paper I considered the CR diffusion zone as an unbounded 
three-dimensional space which is actually too far from the real geometry of the Galaxy. The present work is a continuation of the earlier 
work, but considers the diffusion region as a flat cylindrical disc having both radial and the vertical boundaries. 

In the present study, the propagation of CRs is assumed to follow the same diffusion equation given in Paper I. The solution will be 
applied to local SNRs and the results will be compared to those obtained in Paper I for the burst-like model of particle injection. 



2 CR SPECTRUM FROM A POINT SOURCE 

In the diffusion model, neglecting convection, energy losses and particle losses due to nuclear interactions, the propagation of CR protons in 
the Galaxy is given by the equation 

V ■ (DVN) + Q = ^ (1) 
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where N(r, E, t) is the differential number density, E is the proton kinetic energy, D(E) oc E a with a — constant (positive) is the diffusion 
coefficient which is assumed to be spatially uniform in the Galaxy and Q(r, E, t) is the proton production rate. 

The CR propagation region is assumed to be a cylindrical box bounded in both the radial and vertical directions, and our calculation 
takes into acount the exact location of the sources with respect to the Earth. Inspite of the fact that the actual spatial distribution of observed 
SNRs extent as far as ~ 800pc from the Galactic plane (Stupar et al. 2007), most of the CR propagation studies assume the sources to be 
uniformly distributed in a thin disc of half-thickness ~ (150 — 200)pc. Such an approximation is valid in the study of global properties of 
Galactic CRs since majority of the sources are confined within ~ 200pc from the plane. But, in studies like the present one where the effects 
of nearby discrete sources are discussed, the actual position of the sources should be considered since, for example, for the same source 
distance n(xi,yi, Zi) we expect to see different CR fluxes at different source heights due to the presence of the vertical halo boundary. Our 
calculation will also assume that the Sun is located on the Galactic plane since our Solar system is only ~ 15pc away from the plane (Cohen 
1995). 

The Green's function G(r, r', t, t') of Eq. (1), i.e. the solution for a 5-function source term Q(r, t) = S(r — r')8(t — t') can be found 
so that the general solution can be obtained as 



/OO ft 
dr I dt'G(r,r',t,t')Q(r',E',t') 
-oo J —oc 

Since the CR particles are assumed to be liberated at time t = t\ the equation for G(r, r', t, t') at t > t' becomes simply 



v.( D v G) = f 



(2) 



(3) 



Eq. (3) is solved using the proper boundary conditions and the continuity equations. While solving, we consider the origin to be located at 
(xi, yi,Zi) from the Galactic center. Note that later on this point will represent the actual position of the source with respect to the observer. 
Then, the CR density at a point (x > 0, y > 0, z > 0) due to a point source [which is positioned at (xi,y,, Zi) from the Galactic center] 
with age t, is obtained using Eq. (2) as 
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where R and 77 represent the radial and the vertical boundaries of the Galaxy respectively. The solution at (x < 0, y < 0, z < 0) is obtained 
by just replacing (xj, zi) with {—Xi, —yi, —Zi) in Eq. (4). The proton flux can be calculated using 1(E) ^ (c/4n)N(E), where c is the 
velocity of light and the source spectrum q(E) is taken as 



q(E) = k(E 2 + 2Em p y (T+1)/2 (E + m p ) 



(5) 



in which rn p is the proton mass energy and k is the normalization constant. The source spectral index F is chosen such that F + a = 2.73, 
the observed proton spectral index (Haino et al. 2004). 

For very large radial boundary (R = oo), the solution of Eq. (1) at z > can be written as 
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Fig. 1 compares the proton flux at the Galactic Center given by Eq. (4) with that of Eq. (6) for an SNR-like source located at (1, 0, 0)kpc 
away from the Center with an age t = 2 x 10 4 yr. The results of Eq. (4) at R = (2, 3, 10)kpc are shown by the thin solid, dashed and dotted 
lines respectively. The thick solid lines represent the unbounded solution given by Eq. (6) (i.e. the solution for R — oo). The calculations are 
done at H = lkpc and at H = 3kpc assuming to = 0, represented by the left- and right-hand figures respectively. The diffusion coefficient 
is taken as D(E) = 2 x 10 28 (E /5GeV)°- 6 cm 2 s" 1 for E > 5GeV, where E is in GeV (Engelmann et al. 1990) and the injected protons 
are assumed to carry 10 percent of the total explosion energy of ~ 10 51 ergs. The figures clearly show that, for sources near to the observer, 
the solution of Eq. (4) can be very well approximated by the much simpler unbounded solution for any value of H if 7? > 2>kpc. For example, 
the results at 7? = lOfcpc exactly coincide with the 7? = oo lines. Therefore, considering the fact that our solar system is positioned at a 
distance of ~ 8.5fcpc from the Galactic center and that the Galactic radius extends as far as ~ 20fcpc, the effect of the radial boundary 7? 
on the observed CRs should be negligible at least for those sources that can give appreciable density fluctuations at the Earth, i.e. for those 
sources located within ~ 1.5kpc from the Earth (see Thoudam 2006a). In the following sections where we study the effect of nearby SNRs 
on the observed CRs, we will therefore adopt the simpler Eq. (6) instead of the complicated Eq. (4). 
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Figure 1. CR proton flux at the Galactic center due to an SNR-like source with age t = 2 X 10 4 yr located at (1, 0, 0)kpc away from the center. The thin 
solid line, dashed and dotted lines are the fluxes calculated using Eq. (4) for R = (2, 3, 10)fcpc respectively. The thick solid line represents the flux calculated 
for the boundaryless case (R = oo) using Eq. (6). The calculation assumes to = and D(E) oc E 6 . Left : For H = Ikpc. Right : H = 3kpc. From the 
figures, it can be seen that for any value of H Eq. (4) can be well approximated by the much simpler boundaryless solution if R > 3fcpc. 



3 CR ANISOTROPY 

Knowing the CR density at a point (x, y, z) away from a source of age t, the single source anisotropy amplitude in the diffusion approximation 
can be calculated using (Mao & Shen 1972) 

* = ^"^i (7) 

where Ni is given by Eq. (6) for a point source i located at (xi, yi,Zi) from the Earth. The total anisotropy parameter at the Earth due to a 
number of nearby discrete sources in the presence of an isotropic CR background is given by (Paper I) 

^ ^ li&iVi .Tim 

8 = -* - T (8) 

It 

where the summation is over the nearby discrete sources, r, denotes the direction of the source i giving a flux U and n m denotes the direction 
of maximum intensity. It = 1.37(E/GeV)~ 2 ' 73 cm~ 2 s~ 1 sr~ 1 GeV~ 1 represents the total observed flux of CR protons above ~ WGeV 
(Haino et al. 2004). The phase of the anisotropy is taken as the direction of maximum intensity. Therefore, the anisotropy S as well as the 
phase at an energy E depends on the age and distance of the nearby sources, and may be determined by different sources at different energy 
intervals. However, in the case of a single source dominance, the total anisotropy 8 is given by 8 = (I m / lT)8 m , where m denotes the source 
giving the maximum flux at the Earth. 



4 COMPARISON WITH THE RESULTS OF PAPER I 

In this section, we will try to investigate whether the presence of a halo boundary can affect the anisotropy at the Earth due to nearby sources. 
For that, we consider the 13 known SNRs located within 1.5fepc from the Earth as listed in Table 1 of Paper I. The total anisotropy due to 
these SNRs is calculated using Eq. (8) for different H values at different to's. 

Fig. 2 shows the comparison of the anisotropies calculated in the present work with those obtained in Paper I for the burst-like particle 
injection model. The data points are taken from the compilation of various experiments given in Erlykin & Wolfendale (EW) 2006. Figs 2(a), 
(b), (c) and (d) are the results obtained for to = (0, 10 3 , 10 4 & 5 x 10 4 )j/r, respectively. The thin solid lines represent the results of Eq. (8) 
for H = 0.5kpc, the dashed lines are for H = lkpc, the dotted lines are for H = 2kpc and the dot-dashed lines are for H — 3kpc. The 
thick solid lines are the results of Paper I which were obtained assuming H = oo [Eq. 1 1 of Paper I]. In Figs 2(a)— (c), different sources 
determine the anisotropy at different energy ranges. These are marked by the source names along the lines. It can be seen that the results for 
H — 0.5kpc show a noticeable deviation from the H = oo lines, while those for H = lkpc show a very slight deviation. The results for 
other higher H- values almost overlap with the H = oo lines and are not easily visible in the figures. This shows that, for the particle release 
time of to < W 4 yr, the halo height effect on the local SNR contribution to the observed CR anisotropy is almost negligible if H > lkpc. 
However, the situation is somewhat different in Fig. 2(d) where the calculations are performed at to — 5 x 10 4 yr. Note that this value of 
particle injection time is that at which the model calculated anisotropy values are close to the observed data (see the results of Paper I). The 
anisotropy here is determined solely by the Monogem SNR in the whole energy range considered here, and only those results for H < 2kpc 
show considerable variation from the H — oo line. The results for H > 2kpc show a negligible deviation. Combining all the results of Fig. 
2, we can finally conclude that the effect of the halo boundary of our Galaxy on the local SNR contribution to the observed CR anisotropy 
is negligible as long as the boundary is greater than ~ 2kpc. In the next section, we will combine this result along with the halo heights 
obtained by several authors to discuss the importance of H in the anisotropy study due to local sources. 
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Figure 2. CR anisotropy at the Earth due to nearby known SNRs assuming burst-like injection model. Figs, (a), (b), (c) & (d) are the results obtained for 
to = (0, 10 3 , 10 4 & 5 X 10 4 )yr respectively. Thin solid lines represent the results of Eq. (8) for H = O.bkpc, the thin dashed lines for H = lkpc, dotted 
lines for H = 2kpc and the dot-dashed lines for H = 3kpc. The thick solid lines are the results of Paper I i.e. for H = oo. Figs, (a), (b) & (c) show that 
different sources determine the anisotropy at different energy intervals. These are marked by the source names along the lines. In Fig. (d), the Monogem SNR 
solely determine the anisotropy in the whole energy range. The thick dashed line is the best fit result, in the case of a single source dominance, calculated 
assuming infinite boundaries. Data points are taken from the compilation of various results given in Erlykin & Wolfendale 2006. 



5 DISCUSSIONS AND CONCLUSIONS 

The effect of the nearby CR sources is considered as one of the important effects that can give rise to the observed CR anisotropy at the 
Earth. However, the calculation of CR fluxes from any type of source in the Galaxy essentially requires the use of the proper geometry of 
the Galaxy as well as the actual position of the source with respect to the observer. Since our Galaxy has a cylindrical geometry with the 
radius much larger than the height, the radial boundary is found to have a negligible effect on the CR density and hence the geometry can be 
approximated by an infinite radius with a finite vertical height. Furthermore, this study has found that the effect of the vertical halo boundary 
on the local SNR contribution to the CR anisotropy is negligible if H > 2kpc. 

Fig. 2 shows the effect of the halo height on the CR anisotropy due to nearby known sources for different particle injection times. 
Among the 13 SNRs considered, only Monogem, Vela, G299.2-2.9, SN185 and Cygnus Loop are found to determine the anisotropy at 
different energy intervals. Also, all of them except SN185 (with r = 0.95kpc) have distances r < 0.5kpc. This shows that only the nearest 
sources mainly determine the anisotropy as expected, and hence this results in a negligible halo height effect for H > 2kpc. It is also worth 
mentioning that the vertical heights of the dominant sources above the Galactic plane are found to be less than ~ 60pc which is much less 
than the halo heights {H ^ 0.5fcpc) considered here. 

The actual value of the halo height of our Galaxy is not exactly known. Its value is generally obtained along with other propagation 
parameters using the observed CR data like the secondary /primary ratios, CR density distribution etc. But, the values obtained from the same 
experimental data are different for different CR propagation models. Webber, Lee & Gupta (1992) had obtained a value of H < 4kpc using 
diffusion-convection model. Lukasiak et al. (1994) had obtained H = 2.8^_ggkpc using the Webber et al. (1992) model without convection. 
Webber & Soutoul (1998) obtained H = (2 — 4)kpc and H = (2 — 3)kpc using the diffusion and Monte Carlo models respectively. 
Other results like those of Freedman et al. (1980) and Ptuskin & Soutoul (1998) obtained H ^ 7.8kpc and H = 4.9^ kpc respectively. A 
completely numerical approach using more realistic physical conditions of the Galaxy determined a value of H > 4kpc for the diffusion- 
convection model and H = (4 — 12)kpc for the re-acceleration model (Strong & Moskalenko 1998). These results are found to be consistent 
with the observations of Galactic radio emission structure at 408 MHz which indicate the presence of a thick radio disk with full equivalent 
width of (2.3 ± Q.2)kpc, (3.6 ± QA)kpc and (6.3 ± 0.7) in the Galactic radial range of (0 - 8)fcpc, (8 - 12)fcpc and (12 - 20)fcpc 
respectively (Beuermann et al. 1985), but such a wide range of values makes the Galactic halo height a very uncertain parameter in CR 
propagation studies. However, since most of the values obtained are found to have H > 2kpc, the conclusion given in the previous section 
suggests that the study of local CRs due to nearby SNRs can be carried out without having much information on H. This is because the 
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effect of the nearest sources (r < 0.5fcpc) dominates over the influence of the other nearby sources and the CR fluxes from these sources 
are almost independent of the halo boundary for H > Ikpc as discussed before. Hence, the study of the effect of local sources on the CR 
anisotropy at the Earth can be done using the much simpler three-dimensional unbounded solution. 

For the infinite boundary case, if a single source dominates the anisotropy in the whole energy range as in Fig. 2(d), the total anisotropy 
follows an energy dependence of the form 5 oc E~ a ^ 2 in the high energy regime (Paper I), which for a = 0.6 goes as 8 cx E~ 3 . Such a 
decrease with energy is in fact observed in the high energy anisotropy data somewhere above E ~ 4 x 10 3 GeV upto around 3 x WGeV. 
Moreover, the increase in anisotropy from ~ (1CP 4 — 10~ 3 ) in the (10 2 — 3 x 10 3 )GeV energy range can also be possibly explained by 
a proper choice of (r, t, to) or rather (r, At = t — to) for the single dominant source. We try to estimate the physical parameters of such a 
source that best fit the data. The best-fitting parameters are found to be r = (0.570 ± 0.023)fcpc and A t = (5.343 ± 0.224) x 10 4 yr, and 
the best-fitting line is shown as the thick dashed line in Fig. 2(d). Thus, for to = the source should have an age of t — A t . However, it is 
possible to obtain a number of (t, to) combinations which equally fit the data, all of them giving the same value of r and A t . Therefore, the 
present study only gives an estimate of the distance to the single dominant source; it does not give any precise information on the age and the 
particle release time of the source. It should be noted that it is not the individual (t, to) values that determine the contribution of the source, 
but the propagation time A t of the particles after their release from the source. We can determine the best-fitting t- value only if we know to, 
but the value of to is not exactly known. It may even be that to is an energy-dependent parameter, i.e. particles with different energies emitted 
at different times. Studies based on diffusive shock acceleration in SNRs have shown that the highest energy particles start leaving the source 
region already at the beginning of the Sedov phase (Berezhko et al. 1996), but the major fraction of accelerated CRs remain confined for 
almost around 10 5 j/r for an interstellar medium (ISM) hydrogen atom density of nu = 1cm -3 . This implies that for the local ISM which 
has nn ~ \cmT i (see e.g. Thoudam 2006b and references therein), if a single source determines the whole anisotropy, the source should 
have a characteristic age of ~ 1.5 x 10 5 yr. Unfortunately, there is no nearby known SNR with such an age located at r ~ 0.57fcpc. However, 
it is quite possible that the single dominant source may be an undetected old SNR. In fact, studies assuming adiabatic phase in SNR evolution 
have shown that the surface brightness of an SNR of age ~ 10 J yrs lies below the detection limit of radio telescopes (Leahy & Xinji 1989). 
The present result is further supported by the fact that almost all the nearby sources are quite young with estimated ages less than 10 5 yr (the 
generally accepted particle release time), and they might not have released the CRs into the local ISM. In addition, the possiblity that some 
of the observed features of CRs may be due to undetected nearby sources cannot be simply ignored. 

The single-source explanation of the observed CR properties can also be found in some earlier works (e.g. EW 2000 and references 
therein; EW 2006, etc.), but in a somewhat different context. EW 2000 claimed that the knee in the CR spectrum at E ~ 3PeV can be 
attributed to the presence of a single recent supernova (as yet unidentified) in the local region. On the other hand, EW 2006 tried to explain 
the rise in the anisotropy amplitude as well as the change in its phase near the knee using a single source exploded in the direction from 
the Sun downward of the main CR flux, which are predominantly coming from the inner Galaxy. The latter study considered the source 
parameters as similar to those of the Monogem SNR. Although the single source idea has not been readily accepted by the CR community, 
at the same time there is no reason why it should be just neglected. The present study even points out one more observed property of CRs 
that can possibly be explained by the single source model. 
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